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We report the first results of the beam spin asymmetry measured in the reaction ep — > ep7 at a 
beam energy of 4.25 GeV. A large asymmetry with a sin</> modulation is observed, as predicted for 
the interference term of Deeply Virtual Compton Scattering and the Bethe-Heitler process. The 
amplitude of this modulation is a — 0.202 ± 0.028. In leading-order and leading-twist pQCD, the 
a is directly proportional to the imaginary part of the DVCS amplitude. 
PACS numbers: 13.60.Fz, 14.20.Dh, 24.85. +p 



High-energy lepton nucleon scattering experiments re- 
vealed the quark structure of the nucleon and were em- 
ployed to measure the longitudinal quark momentum 
distributions as well as the quark spin distributions. The 
recently developed formalism [Q-Q for the QCD descrip- 
tion of deeply exclusive leptoproduction reactions intro- 
duces Generalized Parton Distributions (GPDs). They 
carry new information about the dynamical degrees of 
freedom inside the nucleon not accessible in inclusive ex- 
periments. It has been shown that in the Bjorken scaling 
regime (high photon virtuality Q 2 , high energy trans- 
fer at fixed Bjorken scaling variable Xb) the scattering 
amplitude for exclusive processes can be factorized into 
a hard-scattering part (exactly calculable in perturba- 
tive QCD) and a nucleon structure part parametrized 
via GPDs. This is called the "handbag approximation" , 
depicted in Fig. la. In addition to the dependence on 
the parton momentum fraction x GPDs depend on two 
more parameters, the skewcdness £, and the momentum 
transfer t to the baryonic system. In the Bjorken regime 
£ — > xb/(2 — xb), and is related to the momentum imbal- 
ance between the struck quark and the quark that is put 
back into the final state baryon. The momentum trans- 
fer t, which can be varied independently in exclusive pro- 
cesses, measures the transverse momentum transfer be- 
tween interacting quarks. The GPDs reduce to the lon- 
gitudinal parton distributions in the limit £ — ► 0, t — ► 0. 

Deeply exclusive experiments probe the nucleon wave 
function at the quark-gluon level, expressed in terms 
of two spin-dependent and two spin-independent GPDs. 
Deeply Virtual Compton Scattering (DVCS), i.e. elec- 
troproduction of photons from nucleons in the Bjorken 
regime is most suitable for studying GPDs at moder- 
ately high energies. While exclusive meson production 
requires high energies and high photon virtualities to 
reach the Bjorken regime, DVCS may access the GPDs 
already at Q 2 as low as 1 (GeV/c) 2 ||||. 

In this letter we present the first observation of the 
fully exclusive DVCS signal in the beam spin asymmetry 
measured in the reaction ep — > epj. The measurement 
was performed using a 4.25 GeV longitudinally polar- 
ized electron beam and the CEBAF Large Acceptance 
Spectrometer (CLAS) @ at the Thomas Jefferson Na- 
tional Accelerator Facility. Observation of the beam-spin 
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FIG. 1. Feynman diagrams for DVCS and Bethe-Heitler 
processes contributing to the amplitude of ep — > ep7 scat- 
tering. 



asymmetry in semi-exclusive production of high energy 
photons has been reported recently by the HERMES col- 
laboration M , and cross section measurements of DVCS 
by the HI collaboration at DESY at very small xb M, 
which is complementary to the valence quark regime of 
this measurement. 

The experiments measures the DVCS contribution 
through its interference with the Bethe-Heitler process. 
In contrast to DVCS, where the photon is emitted from 
the nucleon, BH photons are emitted from the incom- 
ing or scattered electrons (Fig.l). While the BH cross 
section in most of the kinematic region is much larger 
than the DVCS cross section, the interference between 
the two amplitudes boosts the effect of DVCS and pro- 
duces a large cross section difference for electron helic- 
ities aligned or anti-aligned with the electron beam. In 
this difference, the large helicity-independent BH contri- 
bution drops out, and only the helicity-dependent inter- 
ference term remains [fUJl: 



dQ 2 dxBcttd(j) dQ 2 dxsdtd(j) 
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Here a and b are independent of <j), the angle between 
the lepton and hadron planes, and M 1,1 and M ' 1 are 
the helicity amplitudes for transverse and longitudinal 
virtual photons, respectively. The cross section for posi- 



tive (negative) electron beam helicity is denoted with tr + 
(<x — ). This difference measures directly the imaginary 
part of the DVCS amplitude, since the BH amplitude 
can be calculated exactly. In leading-order only M 1 ' 1 
contributes ]lO| |. 

Experimentally it is much simpler to measure an 



asymmetry A = [d u — d 



-)/(d 4 
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rather than a cross section difference. This asymme- 
try gives the relative strength of the relevant amplitudes 
and has a more complex dependence on <fi than the cross 
section difference, due to <p terms in the denominator. 
In the Bjorken limit, however, the <f> dependence of the 
denominator is expected to be small. 

In this analysis, we used CLAS electroproduction data 
taken in March 1999. Scattering of 4.25 GeV longitudi- 
nally polarized electrons from a liquid hydrogen target 
was studied in a wide range of kinematics. The analyzed 
event sample corresponds to an integrated luminosity of 
1.3 fb" 1 . 

The reaction ep — > epj is identified by analyzing 
the missing mass squared distributions in the reaction 
ep -> epX : M 2 = (u + M- E p ) 2 - (q- P p ) 2 , where v 
and q are the virtual photon energy and momentum, E p 
and P p are the energy and the momentum of the recoil 
proton, and M is the mass of the proton. 

The main background to the single photon final state 
is from 7r° production. At large missing energy the miss- 
ing mass resolution of CLAS is not sufficient to separate 
single photons and ir°s event-by-event. The number of 
single photon events was determined using a fitting tech- 
nique that analyzes the line shape of the M 2 distribu- 
tions. As a fit function the sum of two Gaussians and a 
third order polynomial is used: 

F(M 2 ) = iV 7 • e-( M * 2 - M ^' 2 < 

-{Ml - Mlf/2al 



N^a ■ e 

JZPriM 2 x y 

3=0 



(2) 



The Gaussians represent the M 2 distributions for single 
photon and single pion final states. The polynomial is 
used to model the smooth background that arises from 
radiative processes. The parameters M 2 and a 2 are de- 
termined from the fit to the M 2 distribution of selected 
Bethe-Heitler events (Figure g.a). For M 2 and a 2 , the 
results of a similar fit to the M 2 distribution of selected 
events with final state ep — > ep'fj are used (Figure g.b). 
In these events both photons are detected in the CLAS 
electromagnetic calorimcter:they have an invariant mass 
near the ir° mass. In Figured the curves represent the fit 
to a sum of the Gaussian and a third order polynomial. 
The samples of BH and ir° events were selected in the 
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FIG. 2. Missing mass squared distribution of the detected 
(ep) system for a) ep — > epj and b) ep — > epir . In each 
plot the solid line is the fit to the sum of a Gaussian and 
the third order polynomial distribution. The dashed curve 
corresponds to the Gaussian function and the dotted curve 
represents the polynomial function. 



same kinematical range as the main sample for DVCS 
studies. 

The M 2 distributions are analyzed in eleven <f> bins. 
The final fit to the M 2 distributions with the function in 
Eq. is performed separately for each helicity state and 
for the sum. In the fit only two parameters are adjusted: 
the number of single photon (iV 7 ) and the number of ir° 
(Ntto) events. For each cf> bin the shape of the background 
was fixed from a fit to the points in the summed M 2 dis- 
tribution that excluded the M 2 range from -0.07 GeV 2 to 
0.08 GeV 2 . The relative magnitude of the background 
for distributions at positive and negative helicities are 
found from the fit to the same M 2 range as for the sum. 

Figure || shows a typical fit result for the sum of two 
helicities at <f> = 90°. The solid line corresponds to the 
fitted function F(M 2 ), the dashed line represents the 
Gaussian for the missing photon, and the dashed-dotted 
line is the Gaussian for the missing ir°. The dotted line 
represents the polynomial background. In this bin the 
number of photon and pion events are iV" 7 = 4201 ± 109 
and N^o — 2010 ± 101, respectively. The \ 2 P er degree 
of freedom of the fit is 1.8. Similar results have been 
obtained for all kinematical bins. 

The fitted number of single photon events are used to 
calculate the beam spin asymmetry as: 



A = 
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FIG. 3. Missing mass squared distribution for the reaction 
ep — > epX . Events are integrated in the range of <p from 
70° to 110°. The curves are described in the text. 

where P e is the beam polarization, N-y is the num- 
ber of ep — > ep7 events at positive (negative) beam 
helicity. The average beam polarization, P e = 70%, was 
measured using M0ller scattering. 

In Figure we present our main result, the </> depen- 
dence of A. Data in each 4> bin are integrated in the range 
of Q 2 from 1.00 (GeV/c) 2 to 1.75 (GeV/c) 2 and -t from 
0.1 (GeV/c) 2 to 0.3 (GeV/c) 2 . The error bars shown are 
statistical. Most of the systematic uncertainties related 
to the experiment do not contribute to A. Only the error 
in the measurement of beam polarization, ±1.65%, re- 
mains. There is also a systematic error in the calculation 
of N.y due to the determination of the mean and stan- 
dard deviation of the Gaussian functions for the photon 
and pion missing mass squared distributions, and also 
due to the fit procedure to the M 2 distributions. These 
errors are defined as a deviation of A from its central 
value, when the mean and standard deviation of Gaus- 
sians are shifted within their errors, and when different 
fit techniques were used (see Q for details). 

The data points are fitted with the function A (0) = 
a sin + j3 sin 2<\>. The fitted parameters are a = 0.202 ± 
0.028 s * a *±0.013^ s and/3 = -0.024±0.021 sta *±0.009 S!/s . 
In the Bjorken regime j3 should vanish, leaving only the 
contribution from transverse photons (see e.g. Ref. ]10[). 
In Figure (J the dark shaded region corresponds to the 
range of the fitted function within the statistical uncer- 
tainties of a and (3. The light shaded region includes 
systematic uncertainties on these parameters, estimated 
using the method described above. 

The resulting asymmetry is in a good agreement with 
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FIG. 4. <j> dependence of the beam spin asymmetry A. The 
dark shaded region is the range of the fitted function A((j)) 
defined by the statistical errors of parameters a and /3, the 
light shaded region includes systematic uncertainties added 
linearly to the statistical uncertainties. The curves are model 
calculations according to Refs. [6,11] and are discussed in the 
text. 

a sine/) modulation. Curves in Figure show the re- 
sults of theoretical calculations from Refs. g|ll],[l3|] at 
fixed values of Q 2 = 1.25 (GeV/c) 2 , x B = 0.19, and 
—t = 0.19 (GeV/c) 2 . The limited experimental informa- 
tion does not allow to unambiguously extract GPDs from 
the measurement. Description of GPDs that model the 
£ and t dependencies are therefore used to predict ob- 
servables accessible in experiments. The dashed curve is 
a calculation at leading-twist |5j and no £ dependence in 
the evaluation of GPDs, the dotted curve is leading- twist 
with £ dependence |5| , and the solid curve includes twist- 
3 H,[15| effects. All three calculations include the D- 
term in the parameterization of the GPDs |1(| , which is 
related to double pion contributions. For a more detailed 
description of the model assumptions we refer to a recent 
review [O] . We have estimated that the model asymme- 
tries would be reduced by about 7% if they are averaged 
over the experimental acceptances, bringing them some- 
what closer to the measured data points. 

Although the experimental results are close to the 
lower range of theoretical predictions, none of the cal- 
culations is in agreement with our data. This could be 
due to several factors which need to be studied in fu- 
ture research. Firstly, parameterizations of the £ depen- 
dence of the GPDs were modeled in only a few different 



ways. Secondly, higher order perturbative QCD correc- 
tions have not been included. They tend to reduce the 
measured asymmetry p7[ . Thirdly it is possible that 
the average Q 2 and W of our measurements are not 
sufficiently high to fully reach the Bjorken regime, so 
that some non-perturbative corrections may be needed. 
This can be tested in a future experiment by measuring 
DVCS at the same xb and t but at higher Q 2 and W. 
Such experiments have been approved and are currently 
in preparation |13|h| . 

In conclusion, we have presented the first measure- 
ment of the beam spin asymmetry in exclusive electro- 
production of real photons in the deeply inelastic regime. 
We see a clear asymmetry as expected from the interfer- 
ence of the DVCS and BH processes (EqJ|) . Our results 
agree in sign, and are not far in magnitude to predic- 
tions based on available models of GPD parameteriza- 
tions. This strongly supports expectations that DVCS 
will allow access to GPDs at relatively low energies and 
momentum transfers. This opens up a new avenue for 
the study of nucleon structure, which is inaccessible in 
inclusive scattering experiments. Further measurements 
at higher beam energy are planned, which will allow sig- 
nificant expansion of the Q 2 and xb range covered in 
these studies. The high luminosity available for these 
measurements will make it possible to map out details 
of the Q 2 , xb, and t dependences of GPDs. 
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